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Abstract. We study the sensitivity of multi ton-scale time projection chambers using a 
liquid xenon target, e.g., the proposed DARWIN instrument, to spin-independent and spin- 
dependent WIMP-nucleon scattering interactions. Taking into account realistic backgrounds 
from the detector itself as well as from neutrinos, we examine the impact of exposure, energy 
threshold, background rejection efficiency and energy resolution on the dark matter sensitivity. 
With an exposure of 2001 x y and assuming detector parameters which have been already 
demonstrated experimentally, spin-independent cross sections as low as 2.5 x 10 -49 cm 2 can be 
probed for WIMP masses around 40 GeV/c 2 . Additional improvements in terms of background 
rejection and exposure will further increase the sensitivity, while the ultimate WIMP science 
reach will be limited by neutrinos scattering coherently off the xenon nuclei. 
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1 Introduction 

The nature of dark matter, contributing about 27% to the matter and energy content of our 
Universe [1], is one of the outstanding open questions in physics. A promising dark matter 
candidate is the weakly interacting massive particle (WIMP) [2], which is searched for in 
various experiments world-wide. These projects aim at its direct detection by measuring the 
nuclear recoil signals left by WIMPs scattering in sensitive detectors operated underground [3]. 
Dual-phase time projection chambers (TPCs) filled with liquid xenon (LXe) [4-6] measure the 
light (SI) and the charge signal (S2) generated by a particle interaction. A drift held across the 
TPC removes the ionization charges from the interaction site for detection. The combination 
of both channels allows for the reconstruction of the event vertex and its multiplicity, and is 
used to distinguish signal-like nuclear recoil (NR) events, induced by WIMPs, neutrons and 
neutrinos, from background-like electronic recoils (ER), generated by j3- and 7 -radioactivity 
and by neutrinos interacting with atomic electrons. The high density of LXe provides efficient 
self-shielding, allowing us to define a low-background, central detector region (hducialization). 
Dual-phase TPCs are leading the held in terms of sensitivity since several years [7—11]. No 
dark matter signal has been observed yet [12], excluding spin-independent WIMP-nucleon 
scattering cross sections above 7.6 x 10 _ 46 cm 2 at a WIMP mass of 33GeV/c 2 [10]. 

The next generation of LXe experiments is currently in the design phase or under con¬ 
struction. With target masses beyond the ton-scale, they aim at being sensitive to cross 
sections down to a few xl0~ 47 cm 2 (XENON1T [13], with exposures around 2t x y) or even 
to a few xl0 _ 48 cm 2 (XENONnT [14], LZ [15], 20txy). This poses the question which 
sensitivities could be achieved by the next-to-next generation of instruments, which aim for 
multi-ton target masses and are for example studied within the DARWIN project [16-18]. 
Such detectors will be eventually limited by irreducible background events from coherent 
neutrino-nucleus scattering (CNNS), which are a priori indistinguishable from WIMPs. At 
the same time, all other background sources must be kept under control in order to reach 
the best possible sensitivity. The low-background environment will open up a number of 
additional physics channels, among them the precise measurement of the low energy solar 
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pp-neutrino flux [19, 20], supernova neutrinos [21], axions and axion-like particles [22, 23], as 
well as rare nuclear processes such as the neutrinoless double beta decays of 136 Xe [19] and 
134 Xe, 126 Xe and 124 Xe [24]. However, the search for WIMP dark matter remains the prime 
science case for a DARWIN-like detector. It will be probed via spin-independent and spin- 
dependent WIMP-nucleon interactions, as nat Xe contains about 50% isotopes with non-zero 
nuclear spin, and via inelastic processes [25]. 

In this work we present a detailed study on the sensitivity reach of a DARWIN-type 
multi-ton LXe detector. We mainly focus on spin-independent interactions and take into 
account backgrounds, background rejection efficiency, exposure, thresholds and energy reso¬ 
lution. To this end, we determine the sensitivity, which we define as the average 90% confi¬ 
dence level (CL) exclusion limit, for various combinations of these parameters by simulating 
the outcome of a large number of trial experiments. The analysis is performed by modeling 
all expected signals and backgrounds within a single framework. This allows us to formulate 
a number of requirements which must be considered in the design of such an experiment in 
order to explore cross sections as low as a few xl 0 ~ 49 cm 2 . 

The article is structured as follows: the SI and S2 signal generation in the simulation 
is described in Section 2, while Section 3 explains how we use these signals to construct 
different energy scales, based on the combined S1+S2 signal or the SI signal alone, with their 
impact on the energy resolution. Section 4 describes the method to evaluate the sensitivity 
to spin-independent WIMP-nucleon interactions, based on repeated toy-experiments with 
random background realizations. The background conditions are defined in Section 5 and 
the required ER rejection level is discussed in Section 6 . The main results of this study, the 
dependence of the WIMP sensitivity on various parameters, are presented in Section 7. The 
article closes with a summary and a discussion of our main results, including an extension to 
spin-dependent interactions. 

2 Light and Charge Signal Generation 

In order to simulate the detection process of signal and background events in a realistic 
fashion, we employ a signal generation model following Ref. [26]. Based on the energy and 
the type of interaction, ERs measured in keVee (electronic recoil equivalent energy) or NRs 
measured in keVnr (nuclear recoil equivalent), the mean number of generated quanta n 7 
(scintillation photons) and n e (ionization electrons) is calculated for a drift field of 500 V/cm, 
a typical number for dual-phase TPCs. For ERs, charge and light yields are taken from NEST 
v0.98 [26], while for NRs, Q y and £ e g- are taken from XENON100, using Refs. [27] and [28], 
respectively. Q y is only given down to 3 keVnr and we extrapolated to lower energies by 
fixing it to the value at 3 keVnr. The £ e ffparametrization assumes £ e g- = 0 for E < 1 keVnr, 
with the consequence that no light is produced for energy depositions below this energy. The 
respective event is therefore discarded, as it would not be observed in the TPC. The total 
number of quanta n = n 7 + n e is convoluted by a Gaussian distribution with a width of 
<7 = y/F x n, assuming a Fano factor F = 0.03 [29], in order to get the generated number of 
quanta q. The number of scintillation photons, eventually leading to the light signal (SI), is 
randomly drawn from a Binomial distribution 

n 9 = Binomial(g,p = ~zr)- (2-1) 

' n 

The observable in the detector is the number of photoelectrons (PE) 

n° s i = n 7 x e 
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with the photon detection efficiency 
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L y is the average light yield in PE/keVee at 122 keV and zero held and Y = 63.4 is the number 
of photons/keV generated in LXe at zero held [26]. Our benchmark value is L y = 8PE/keV, 
about 2x larger than achieved in XENONIOO [60] and similar to the LUX value [61]. The 
number of electrons from the charge signal (S2) is consequently 


n° S 2 = Q ~ n ° r ( 2 - 4 ) 

n ° s j and n^ 2 are anti-correlated, and the resolution due to the statistical nature of the signal- 
generation process is properly modeled. The signal quanta observed in the detector, nsi and 
ns 2 , are drawn from a Poisson or a Gaussian distribution, respectively, with the means 
and n^ 2 and resolutions tuned to match existing experimental data. This takes into account 
additional effects from the photoelectron detection process (“single photoelectron resolution”) 
as well as the conversion of ionization electrons to photoelectrons (Gaussian process with a 
mean of 20-30 PE/e _ [30, 31]). Throughout this study, the charge signal (S2) is given in 
electrons, and we do not consider charge losses during the drift of the electrons through the 
LXe, which requires an effective gas cleaning system and 100% charge extraction into the gas 
phase (as achieved by XENONIOO [30]). 


3 Energy Scale and Energy Resolution 

In this work we consider two energy scales. The first one relies on the scintillation signal (SI) 
only, which is still the standard in the held and used in most published analyses [8, 9, 11]. The 
second one is based on light and charge (S1+S2), additionally taking into account the anti¬ 
correlation between the two observables, which is induced by the signal generation process: 
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Figure 1. (left) Simulated light (SI) and charge signals (S2) from hypothetical mono-energetic NR 
energy depositions. The same number of events was simulated for all three energies. A light yield 
L y = 8.0 PE/keVee, full charge extraction, and the processes discussed in Section 2 are considered. 
The anti-correlation between both signals is very weak, (right) Resolution of different energy scales 
as derived from the simulation. The “combined energy” scales use the SI and S2 signals and take into 
account the anti-correlation. The data points are well described by functions of the form a o + ai/y/E. 
As the light signal is non-symmetric at low energies, the RMS is a better measure for the resolution 
than the Gaussian width a. 
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ionization electrons which recombine with Xe^ molecules will be lost to the charge signal 
but eventually create scintillation light. This anti-correlation, though well known for high- 
energetic 7 -lines [32], has not yet been observed directly for low-energetic nuclear recoils. 
The event generation model described above allows us to establish the combined energy 
scale by simulating the light and charge signals generated by discrete NR energy depositions 
in the LXe, see Figure 1 (left). The anti-correlation is only very weak and hardly visible 
when the signal resolution is applied correctly. The energy scales are constructed based 
on the mean observed signals obtained for many mono-energetic energy depositions. Using 
the combined energy scale, which additionaly employs an anti-correlation angle, leads to a 
significant improvement in the energy resolution compared to the Sl-only scale. However, the 
improvement compared to a charge-only scale is small, see Figure 1 (right). It is interesting 
to note that a L y increase from 8 PE/keVee to 12PE/keVee does not lead to a significant 
improvement in the resolution of the combined scale. 

We will see in the following that the CNNS background induced by solar 8 B-neutrinos 
increases exponentially at very low NR energies. A good energy resolution is therefore manda¬ 
tory in order to decrease the energy threshold for the WIMP search to values as close as 
possible to the rise. For an Sl-based scale ( 8 PE/keVee), one can achieve 6 keVnr, while a 
threshold of 5 keVnr or lower is only reachable when a combined energy scale is applied. 

4 Evaluation of the spin-independent WIMP-nucleon Sensitivity 

We evaluate the sensitivity of a DARWIN-type multi-ton LXe detector to spin-independent 
WIMP-nucleon scattering cross-sections by averaging the outcome of many simulated pseudo¬ 
experiments. The analysis is performed under realistic assumptions, i.e., the statistical nature 
of the signal generation (energy resolution) and of the observed background conditions is 
taken into account. All energy depositions in the detector, from signal and background, are 
converted into a reconstructed NR energy. For ER background events of ER energy E ee , the 
respective SI and S2 signals are generated according to Section 2. From these observables, 
the reconstructed NR energy is calculated. For NR background events of energy E nr . SI and 
S2 signals are generated as well and also back-converted into reconstructed NR energy. 

For each pseudo-experiment, the (integer) number of background events and their en¬ 
ergies in reconstructed NR energy are determined by independently considering each of the 
background sources X introduced in Section 5 below. For an expected mean number of back¬ 
ground events bx ■ for a given exposure, energy interval and ER rejection level, the actual 
number of events nx is randomly drawn from a Poisson distribution with mean b\■ The 
energies of these events are distributed according to the energy spectra defined in Section 5, 
taking into account the energy scale under study. This means that background events from 
outside the WIMP analysis window, defined in reconstructed NR energy, can enter the region 
of interest (ROI) and events with initial energies in the window can be lost if the reconstructed 
NR energy is not in the ROI. The total number of background events is rib- 

The expected spectra for WIMPs interacting in a LXe target of total mass M are cal¬ 
culated from the differential event rate 

dR po M [°° , da s 

-rPT- = - / 'yf(v)—- dv. (4.1) 

dE nr m Xe m x J Vmin dE nr 

rrixe an d m x are the masses of the xenon nucleus and the WIMP, respectively, and f(v) is the 
normalized WIMP velocity distribution. We take into account the masses and abundances of 
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the different xenon isotopes occurring in natural xenon gas. All velocities are defined in the 
detector’s reference frame, with 


'Vmin 


'E nr m Xe (rn Xe + m x ) 2 


(m.v e m x ) s 


E nr HlXe 1 




(4.2) 


being the minimal velocity required to induce a nuclear recoil E nr . The velocity distribution 
f(v) is truncated at the escape velocity v esc = 544kms _1 [33], the maximum velocity of 
WIMPs bound in the potential well of the galaxy. The canonical value for the local WIMP 
density at the Sun’s position in the Galaxy is po = 0.3GeV W 2 cm” 3 . Because of its large 
de Broglie wavelength, the WIMP interacts coherently with all nucleons in the target nucleus. 
We mainly focus on spin-independent (SI) interactions for which the WIMP-nucleus scattering 
cross section in Eq. (4.1) is given by 


da 

dE nr 


m N 

2uV 


asiF 2 (E nr ). 


(4.3) 


The loss of coherence for heavy WIMP targets such as xenon is accounted for by the finite 
form factor F(E nr ). for which we use Helm’s definition. The cross section reads 
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(4.4) 


where the f p ^ n describe the WIMP couplings to protons and neutrons. The second equality 
assumes f p = f n , leading to an A 2 dependence of the cross section, fi is the WIMP-nucleus 
reduced mass, see Eq. (4.2), and fi n the one of the WIMP-nucleon system. It is used to 
relate the WIMP-nucleus cross section a to the WIMP-nucleon cross section a n , which is the 
relevant parameter for this study and allows the comparison with other target nuclei. 

For a single event, the reconstructed NR energy E' nr is obtained by generating the 
corresponding SI and S2 signals and the subsequent back-conversion using the previously 
defined energy scale. The observed number of events in a simulated toy experiment running 
for a live-time T is obtained by integrating Eq. (4.1) from the threshold energy Et to the 
upper boundary E' high : 



In most parts of this study, we take an energy-independent efficiency e = 0.3 for NRs. As 
dR/dE nr is a steeply falling exponential function, E' high is much less relevant than the energy 
threshold Ei and fixed to 20.5keVnr or 35.0keVnr. 

The WIMP sensitivity is evaluated in two ways: in order to study how the sensitiv¬ 
ity depends on various parameters we perform an analysis based on a fixed WIMP search 
window, which is defined from E[ to 20.5keVnr and by a fixed ER rejection level at 30% 
NR acceptance. For each simulated trial, all background events in the box are considered 
as possible WIMP signals, and the 90% CL exclusion limit is calculated according the maxi¬ 
mum gap method [34], By construction, this method takes into account the different spectral 
shapes of signal and background. The sensitivity is determined by averaging the results of 
1000 individual simulations. The advantage of this approach is that it is straight-forward to 
quantify the number of signal and background events in the region of interest, and allows for 
an easy variation of detector parameters. The upper energy limit of 20.5keVnr was chosen 
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as the signal-to-background ratio starts to deteriorate at higher energies, especially for small 
(~10GeV/c 2 ) and medium WIMP masses (~100GeV/c 2 ), and the gain in sensitivity is small. 

Systematic uncertainties, such as the finite knowledge on the energy scale which should 
be considered in a proper 90% CL exclusion limit, are beyond the scope of this study and 
neglected. It is known that likelihood methods, also taking into account the event distribu¬ 
tion in discrimination space DocS2/Sl, can improve the sensitivity by a factor ~2 [35] and 
additionally include systematic uncertainties in a proper way. For the optimal parameters 
found by the study with the fixed WIMP search window, we calculate the sensitivity using a 
likelihood method. The natural logarithm of the likelihood is defined as 

TLtot 

log £ = b ER + &iVi?) + ^2 log [ Vxfx( E nr,k) 9x( D k) + bER fER{E' nrk ) gER(D k ) 

k =1 

+ b NR f NR {E' nrk ) gNR(D k ) ], (4.6) 

where /i x and b\ are the expected mean values for WIMPs and ER and NR backgrounds. 
fx{E ' nrk ) and gx{D k ) denote the normalized probability distribution functions of event 
source X in energy E' nr and discrimination parameter space D , respectively. The sum is 
evaluated for all for all ntot events generated in each toy experiment. The test statistic is 
defined as 

g=/ - 2 l 0 g (£fer) = 21 og£(// x )- 21 og£(// x ) for g x < g x ^ ^ 

lo for fi x > fi x , 

where fi x is Maximum Likelihood estimator of Eq. (4.6). Its distributions under the signal 
and background hypothesis Hq are used to derive the mean sensitivity. The procedure 
takes into account the CLs mechanism to protect the result against statistical downward 
fluctuations of the background [36]. We generally follow [35, 37] where more details can be 
found. 

5 Backgrounds 

We consider all relevant backgrounds for multi-ton scale LXe dark matter detectors, which 
can be separated into two categories. External backgrounds comprise 7 -rays and neutrons, 
which stem from radioactive decays or interactions outside of the LXe target, for example 
in the detector’s construction materials. Due to the high density of LXe (p ~ 3g/cm 3 ), 
they can be reduced considerably by target fiducialization. This does not work for intrinsic 
backgrounds, which are uniformly distributed in the target region. 

In the following, we detail the background assumptions which enter our analysis. All 
contributions leading to single scatter signatures in the low-energy WIMP search region are 
summarized in Table 1. If the background distribution is not flat, we quote the average values 
for a given energy interval. All numbers are given assuming an infinite energy resolution while 
for the sensitivity study presented below, the resolution is taken into account. 

7 -rays from materials The 7 -background produced in the laboratory is virtually zero 
behind a shield of several meters of water [14]. The remaining background stems from ra¬ 
dioactive contamination (mainly from the 238 U, 235 U and 232 Th chains, 40 K, 60 Co) in the 
cryostat and detector materials. It can be reduced by target fiducialization and discrimina¬ 
tion using the charge-to-light (S2/S1) ratio. Based on earlier studies [19], we assume a flat 
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Source 

Rate 

[events/(t-y-keV xx)[ 

Spectrum 

Comment 

7 -rays materials 

0.054 

flat 

assumptions as discussed in text 

neutrons* 

3.8xl0 -5 

exp. decrease 

average of [5.0-20.5] keVnr interval 

intrinsic 85 Kr 

1.44 

flat 

assume 0.1 ppt of nat Kr 

intrinsic 222 Rn 

0.35 

flat 

assume 0.1 ^Bq/kg of 222 Rn 

2zv/3/3 of 136 Xe 

0.73 

linear rise 

average of [2-10] keVee interval 

pp- and 'Be u 

3.25 

flat 

details see [19] 

CNNS* 

0.0022 

real 

average of [4.0-20.5] keVnr interval 


Table 1. Background contributions considered in the sensitivity analysis. All values given here 
are before ER discrimination or a finite NR acceptance and do not yet take into account an energy 
resolution. The ER background rates are given per keVee (electron recoil equivalent). The NR 
backgrounds, marked by the asterisk, are given per keVnr (nuclear recoil equivalent). If the background 
is not flat, an average value over a finite interval is quoted. The energy resolution is taken into account 
for the sensitivity study. This is of particular importance in case of the steeply rising CNNS spectrum. 


background spectrum and a rate of 0.054 events/(t-y-keVee) before discrimination, which is 
reached behind a LXe layer of ~13cm. a- and /3-backgrounds from detector materials are 
irrelevant even for small LXe TPCs as they do not penetrate into the fiducial target region. 

Neutrons As WIMPs are expected to generate nuclear recoils in the detector, background 
from neutrons, producing the same signature, is dangerous as the S2/S 1-background rejection 
cannot be used. However, one can rely on target Lducialization (even though neutrons have 
a considerably longer mean free path than 7 -rays of the same energy, and have also higher 
energies) and on rejection based on the event multiplicity. Background neutrons are either 
of cosmogenic origin, i.e., produced in muon induced electromagnetic and hadronic showers, 
or radiogenic, i.e., generated in (a,n) or spontaneous fission of heavy isotopes. A shielding 
of several meters of water, operated as a Cerenkov muon veto, decreases the muon-induced 
neutron background to very low levels [14], and radiogenic neutrons from the laboratory are 
efficiently reduced by many orders of magnitude. Both backgrounds are therefore considered 
negligible in this study. 

More relevant are radiogenic neutrons emitted from the detector materials itself. The 
light PTFE, usually used as insulator and light reflector [38], is a known source of ( a,n) 
neutrons, where the a-particle is from 238 U, 235 U and 232 Th-chain contaminations. An¬ 
other critical component are the photosensors which are made from a variety of materials, 
sometimes without low-background alternatives [39]. Here we assume a NR background 
from radiogenic neutrons which is dominated by the photosensors and contributes a rate of 
3.8 X 10 ~ 5 events/(t-y-keVnr) (average value in a 5.0-20.5keVnr interval). This is compara¬ 
ble to the ER background from two-neutrino double-beta decays (2i/j3j3) when taking into 
account realistic ER rejection and NR acceptance levels. This neutron background rate is 
motivated by a Monte Carlo simulation performed for an an idealized DARWIN-type detec¬ 
tor with a 181 LXe target, where the relevant components are assumed to be made from 
OFHC (cryostat vessels, electrodes, diving bell), PTFE (TPC) and quartz (photosensors). 
We have optimized the detector model presented in Ref. [19] for a low neutron background by 
reducing the TPC PTFE reflectors to a thickness of 5 mm, and by assuming cleaner materials 
(PTFE factor 4, e.g., entry 43 in [40]; copper factor 2, e.g., the “cryostat copper” in [41], and 
photosensors factor 5). These would have to be identified in dedicated material screening 
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and selection campaigns. The contributions of the photosensors could possibly be reduced by 
using improved low-background sapphire instead of quartz [39] or by employing novel sensors 
such as gaseous photomultipliers (GPMs) [42], 

The neutron production rates and energy spectra were calculated using the SOURCES- 
4A code [43], based on the assumed radioactive contamination of the materials. Details 
on the calculations can be found in [44], For each component (the two cryostats and the 
two photosensor arrays were simulated independently), 4 X 10' neutrons were simulated in 
Geant4 [45], which leads to a negligible statistical uncertainty. The finite efficiency to identify 
multiple scatter signatures increases the true single scatter background by a factor 1.85 when 
taking into account a resolution of 3 mm (10 mm) to separate vertices in the z ( xy ) coordinate 
and assuming a realistic lower threshold to identify the second S2 signal in double scatters [44] . 
Figure 2 (left) shows the exponentially falling recoil spectrum, which is reached after a fiducial 
cut of ~16cm from all sides. As the surface-to-volume ratio decreases with increasing detector 
size, the background situation will improve for detectors larger than the simulated one, as the 
background sources scale with the surface. For this study, we do not pick a specific geometry 
but assume that the neutron background scales with exposure, similar to target-intrinsic 
backgrounds. We note that a full study of this background contribution requires the precise 
knowledge of the detector components for the neutron production (contamination, material 
assembly), as well as detector size and geometry, which is not defined at this point. 



Figure 2. (left) Differential energy spectrum of single-scatter NR background from the detector 
materials, assuming 100% NR acceptance and an infinite energy resolution, (right) Differential 
energy spectra of the ER background sources. No energy resolution or ER rejection is applied at this 
point. The dominating contribution is from pp- and 'Be solar neutrinos. 


Intrinsic 8o Kr Xenon does not have long-lived isotopes besides 136 Xe, a 2vf3/3 emitter, and 
can be purified to a high degree from other contaminations. One exception are other noble 
gases, as these are chemically inert (cannot by removed by chemical methods) and mix very 
well into xenon. As a consequence, such backgrounds are uniformly distributed in the target. 
A problematic isotope is the anthropogenic /3-emitter 85 Kr (Ijy 2 = 10.76y), present in natural 
krypton at the 2 X 10~ n g/g level [46]. Krypton can be separated from xenon by cryogenic 
distillation, exploiting its 10 X higher vapor pressure at LXe temperatures [47], and nat Kr- 
concentrations of less than lppt in xenon have already been achieved [48]. At low energies, 
the single scatter ER spectrum from 85 Kr is basically flat [49]. For this study we assume a 
natRr concentration of 0.1 ppt, only a factor of ~5 below the design goal of systems currently 
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under development for the next generation projects [50]. It leads to 1.44events/(t-y-keVee) 
before S2/S1 discrimination [19], see also Figure 2 (right). 

Intrinsic 222 Rn The radioactive noble gas 222 Rn is part of the 238 U chain. It decays with a 
short half-life of 3.8 d, however, it is constantly re-emanated from surfaces due to the presence 
of 226 Ra. Therefore, all materials in contact with xenon anywhere in the detector system must 
be selected for their low Rn-emanation [51]. The DARWIN background needs to be around 
0.1 ^rBq per kg of LXe, about an order to magnitude below the goal of the next generation 
experiments. The resulting single scatter rate before discrimination is 0.35 events/(t-y-keVee) 
with a flat recoil spectrum [19]. While the improved surface-to-volume ratio is generally 
favorable for the reduction of this background in very large detectors, an active removal of 
222 Rn from the LXe target is difficult at these low concentrations and would require further 
R&D. Adsorption on charcoal [52], for example, would need enormous amounts of additional 
xenon and charcoal with a sufficiently low 226 Ra contamination. 

Two-neutrino double-beta decay of 136 Xe We assume that the DARWIN facility is 
initially operated with natural xenon, which contains 8.9% of 136 Xe known to decay via a 
two-neutrino double-beta process (2vf3(3) with a half-life of Tj/ 2 = 2.17 x 10 21 y [53]. Even 
though this is 11 orders of magnitude longer than the age of the Universe, its continuous 
single-scatter spectrum (Q-value is (2458.7 ± 0.6) keV [54]) contributes to the WIMP search 
background. At lowest energies, its ER spectrum can be well approximated by a linear 
function starting from (0,0). The integral rate in an 2-10keVee interval is 5.8events/(t-y) 
before ER rejection [19]. In principle this background can be completely avoided by using a 
target depleted of 136 Xe. 

Low-energy solar neutrinos The vast majority of neutrinos emitted from the Sun are 
neutrinos generated in the pp-fusion process or the subsequent "Be reaction. Due to their 
rather low energies and high abundance, together with the impossibility to reduce their con¬ 
tribution by target purification, hducialization or single-scatter selection, they are the most 
relevant source of ER background for LXe-based dark matter searches beyond the ton-scale. 
At lowest energies, they generate a basically flat ER spectrum and can only be separated from 
WIMP-induced NRs based on their S2/S1 ratio, albeit with a finite efficiency. The combined 
rate of pp- and 'Be neutrinos, which are indistinguishable in a DARWIN-like detector, is 
3.25events/(t-y-keVee) before discrimination [19]. Extending the energy range beyond the 
WIMP search region to ~30keVee allows the precise measurement of the low-energy solar 
neutrino flux with DARWIN. 

Coherent neutrino-nucleus scattering Dark matter WIMPs are expected to produce 
single-scatter NR events in the LXe detector. The identical signature is generated by neu¬ 
trinos, which are predicted to coherently scatter off the xenon nucleus. Therefore, it is 
impossible to a priori distinguish the two sources of events. As coherent neutrino-nucleus 
scattering (CNNS) is a standard model process, and neutrinos of the right energies to pro¬ 
duce signals in the WIMP search energy region exist, this background will ultimately limit 
the sensitivity of any WIMP detector [55], if no additional information on the direction of the 
NR track is available [56]. For a LXe target, the relevant neutrinos are mainly 8 B neutrinos 
from the Sun, generating a “wall” of events rising steeply below ~4keVnr, which significantly 
reduces the sensitivity to 5-8 GeV/ c 2 low-mass WIMPs. At higher NR energies, the main 
CNNS background is from atmospheric neutrinos, with smaller contributions from solar neu¬ 
trinos from the helium-proton reaction (hep) and the diffuse supernova neutrino background 
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Figure 3. (left) Differential CNNS spectrum for different energy scales, based on Sl-only or S1+S2. 
The largest improvement in energy resolution can be achieved by using light and charge signal; the 
impact of a light yield increase by 50% is only mild. The spectrum with infinite resolution (dashed) 
is taken from [57]. The steep rise below 4.0keVnr is from solar 8 B neutrinos. We assume that no 
light is produced for E nr < 1 keVnr, leading to the large number of entries in the lowest bin of the 
reconstructed spectra. Using the S1+S2 scale with 8 PE/keVee, we also show the expected recoil spec¬ 
tra for WIMPs of 10GeV/c 2 and 40GeV/c 2 with cross sections of 2 x 10 _ 46 cm 2 and 2 x 10 _ 48 cm 2 , 
respectively, (right) WIMP exclusion limits at 90% CL ignoring all backgrounds besides the one 
from CNNS for exposures of 100 1 x y, 200 1 x y, 5001 x y and 1200 1 x y. For the likelihood analysis 
we use a S1+S2 combined energy scale with a light yield of 8 PE/keVee at 122keVee, a 5-35 keVnr 
energy interval and assume an unrealistic NR acceptance of 100%. (The ER rejection will lower this 
value.) The result for 2001 x y is shown with its la and 2a intervals together with the published 
limits from XENON100 [9] and LUX [61]. For comparison the “WIMP discovery limit” (red dashed) 
and the “1 event line” (red dotted) of [55] are shown as well. 


(DSNB) [57]. The background from CNNS strongly depends on the threshold, see Figure 3 
(left), which renders the energy scale (resolution) used for the analysis crucial. At higher NR 
energies, the rate is fairly flat around 1-2 xlCU 3 events/(t-y-keVnr). 

We note that the “WIMP discovery limit”, as introduced in Ref. [55] for every WIMP 
mass, corresponds to the cross section at which a WIMP can be detected at 3 a given a 
background of 500 CNNS events above a LXe threshold of 4 keVnr (no energy resolution 
applied). A utopic LXe exposure of 53001 X y is required to reach this number of background 
events in a 4-35 keVnr interval, even assuming 100% NR acceptance. Being somewhat more 
realistic, we calculate the average 90% CL exclusion limit for exposures from lOOtxy to 
12001 x y, considering only the CNNS background assuming 100% NR acceptance. We use 
the combined energy scale of Section 3 with a light yield of 8.0 PE/keVee and an energy range 
of 5-35 keVnr. The result is shown in Figure 3 (right) and compared to the “WIMP discovery 
limit” and the “1 event line” of [55]. 

6 Electronic Recoil Rejection 

In order to reach sensitivities which are limited by the “ultimate” NR background from CNNS, 
the ER backgrounds from intrinsic radioactive impurities, two-neutrino double-beta decay and 
especially solar neutrinos interacting with atomic electrons have to be reduced significantly. 
The lifetimes of the xenon singlet and triplet excimer states, which produce the scintillation 
light, are not very different, rendering pulse-shape discrimination inefficient [59]. Therefore, 
LXe-based dual-phase TPCs rely on signal-background discrimination based on the different 
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Light Yield [PE/keVee] 


Figure 4. Illustration of the impact of the zero- 
field light yield L y at 122keVee on the ER rejec¬ 
tion, assuming a constant NR acceptance of 30%. 
Using the signal model introduced in Section 2, 
the fraction of ER events leaking into a low en¬ 
ergy WIMP search region is determined in a sim¬ 
ulation, and is quoted relative to the leakage for 
a moderate light yield of 4.0 PE/keVee. All other 
parameters which might have an impact on the 
ER rejection were kept constant. The fluctuations 
in the leakage fraction are purely statistical. 


charge-to-light ratio (S2/S1) for ERs and NRs, caused by the different energy loss mecha¬ 
nisms [58]: at a given energy deposition, ERs exhibit a larger S2 pulse than NRs. Rejection 
levels of 99.5% [60], 99.6% [10] and 99.987% [8] at 50% NR acceptance have been already 
achieved in dark matter detectors, where the rejection appears to improve towards lower re¬ 
coil energies. In this study, we assume flat, average values for rejection and acceptance. We 
also ignore possible pathological signatures which could lead to so-called “anomalous” leakage 
into the WIMP signal region. Such signatures could come from incomplete light or charge 
collection or from accidental coincidences of causally non-connected pulses. 

The mean values of the S2/S1 vs. energy distributions for ERs and NRs are determined 
by the signal generation processes, which are affected by the drift field in the TPC: a lower 
drift field leads to more charge recombination and hence a smaller S2 but a larger SI signal. 
The width of the distributions, crucial to achieve high ER rejection at high NR acceptance, 
is determined by the intrinsic fluctuations in the generation of the initial quanta and the 
statistical fluctuations in the detection processes. The latter can be mainly addressed by the 
design of a multi-ton scale detector, for example by maximizing the light and charge yields. 
By using the signal generation model introduced in Section 2, we illustrate the impact of the 
light yield on the ER rejection in Figure 4. We estimate the fraction of ER events leaking 
into a low energy WIMP search region, which is defined by a fixed 30% NR acceptance. A 
moderate light yield of 4.0 PE/keVee, as realized by XENON100 [60], is used as a reference 
with its leakage fraction normalized to unity. An increase to L y = 8.0 PE/keVee reduces the 
leakage fraction by almost a factor 4, the higher L y = 12.0 PE/keVee by more than a factor 7. 
All other parameters which might affect background discrimination are kept constant here. 

In order to determine the “typical” ER rejection level required for a multi-ton scale LXe 
dark matter detector (and the corresponding NR acceptance), we compute the numbers at 
which the ER background from solar neutrinos is 20% of the one from CNNS NRs in our low- 
energy window. An energy scale has to be selected in order to evaluate the ER background 
in the NR energy region. Due to the steeply rising CNNS spectrum, the rejection level 
depends crucially on the low energy threshold energy E[ (see also Section 4); the upper 
limit is again fixed to E' high = 20.5keVnr. For a variable NR acceptance, Figure 5 (left) shows 
the required ER rejection level for five different combinations of energy scales, light yields, 
and threshold energies. The resulting typical numbers are around 99.98% ER rejection and 
30% NR acceptance (dashed lines), which we will further use as reference values. We note 
that the achievable WIMP sensitivity depends on the absolute size of the background, which 
varies by more than a factor 2 for the five cases in Figure 5, when evaluated at a fixed NR 
acceptance of 30%. Figure 5 (right) separately shows the mean expected number of events 
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Figure 5. (left) ER rejection level required to achieve an ER background from solar pp- and 'Be 
neutrinos which is 20% of the NR contribution from CNNS, as a function of NR acceptance. The 
comparison in a fixed low energy window (A' — 20.5 keVnr) requires the choice of an energy scale and 
depends on the low energy threshold E' low , due to the steeply rising CNNS spectrum. The figure 
shows five representative cases, the dashed lines indicate the “reference” numbers used in this study, 
(right) Mean number of expected NR background events from CNNS and materials for different NR 
acceptances and mean number of expected ER background events (solar neutrinos, 85 Kr, 222 Rn, 2p/3f3 
and from materials as discussed in Section 5). The solar neutrinos are also shown independently. The 
NR background dominates for ER rejection levels better than 99.9%. We assume an exposure of 
2001 x y, 5.0-20.5 keVnr, and a combined energy scale with 8PE/keVee. 

from all ER (neutrinos, 0.1 ppt of nat Kr, 0.1^Bq/kg of 222 Rn, materials) and NR sources 
(CNNS, materials), where the latter are given for different acceptances. The individual ER 
contribution of solar neutrinos is shown as well. We assume an exposure of 2001 X y, a 5.0- 
20.5 keVnr energy window and a combined energy scale with L y = 8.0PE/keVee. The NR 
background exceeds the one from ERs for rejection levels >99.9%. As the background of an 
ultimate detector should be dominated by CNNS events, the ER rejection must be >99.98% 
at a NR acceptance of 30-50%. 

7 WIMP Sensitivity 

The following experimental parameters were individually varied in a systematic way in order 
to evaluate their impact on the dark matter sensitivity: 

• exposure: target mass M x live-time T, 

• lower energy threshold E', , 

• ER background rejection (at a fixed NR acceptance of 30%), and 

• energy scale: combined (S1+S2) and Sl-only scale using two different light yields. 

The parameters not under study were fixed to the following values, which were found to be 
optimal for a realistic experiment: an exposure of 2001 X y, a lower and an upper threshold 
of 5.0 keVnr and 20.5 keVnr, respectively, and an ER rejection level of 99.98% (2 X 10 -4 ). A 
combined energy scale (CES) with a light yield of L y = 8 PE/keVee was used for the standard 
analysis and all background sources discussed in Section 5 were taken into account. 
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Figure 6. Average sensitivity, defined as the 90% CL exclusion limit on the cross section, at 
~40GeV/c 2 (left column) and at 10GeV/c 2 and 1000 GeV/c 2 (right column). It is plotted vs. (1 st row) 
exposure, (2 nd ) lower threshold E ' low , (3 rd ) ER rejection level and (4 th ) energy scale. All values are 
normalized to the marked reference value, obtained by using “standard” parameters. The parameters 
not under study were fixed to the values indicated in the plots. Open symbols indicate the sensitivity 
considering only the CNNS background (at 30% NR acceptance). See text for more discussion. 
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The results are presented in Figure 6. The four panels in the left column show the average 
sensitivity at the minimum of the sensitivity curve, typically around 40GeV/c 2 , obtained for 
a fixed WIMP search box. The limit at 100GeV/c 2 is only about 60% weaker, rendering this 
choice as typical for medium (weak scale) WIMP masses. The four panels on the right show 
the sensitivity at low (m x = lOGeV/c 2 ) and high (1000GeV/c 2 ) WIMP masses. Fluctuations 
in the curves are of statistical nature. All sensitivities are given relative to the one obtained 
for the standard parameters, marked in the plots by the double ring. 

Exposure We study the WIMP sensitivity for exposures from 1001 x y to 12001 X y (first 
row in Figure 6). The CNNS background becomes relevant above 3001 x y, where it con¬ 
tributes ~3 events to the background in the WIMP search box and the sensitivity starts to 
flatten. The gain in sensitivity is only a factor 1.3 when increasing the exposure from 3001 x y 
to 12001 x y. At m x ~ 40GeV/c 2 , the improvement from the realistic background situation 
to the CNNS-only case (dashed lines, open markers) increases from a factor ~1.05 at 1001 x y 
to ~1.25 at 12001 x y. With an increase from a factor 1.2 to 1.3, the situation is almost iden¬ 
tical at 1000GeV/c 2 . As the background is dominated by CNNS at low recoil energies, the 
difference is less pronounced at 10GeV/c 2 , where the other background sources are basically 
irrelevant, such that the two sensitivities almost agree at 10GeV/c 2 . The gain in sensitivity 
for the CNNS-only case shown here, assuming no ER background but a finite NR efficiency of 
30%, is less pronounced than in Figure 3 (right), where 100% NR was assumed to illustrate 
the most optimistic background scenario. 

Lower energy threshold Changing the lower energy threshold E', has the largest impact 
at small WIMP masses, due to their expected steeply falling recoil spectrum. This is shown 
in the second row of Figure 6, where the sensitivity at 10GeV/c 2 decreases by two orders 
of magnitude when increasing the lower threshold from 3.5keVnr to 9keVnr in a 200txy 
run. Around E[ ~ 5keVnr, the sensitivity starts to level off, and there is no significant 
sensitivity gain by going to lower thresholds, as the irreducible CNNS background effectively 
builds a “wall” of background events, limiting the energy range with an acceptable signal-to- 
background ratio. The sensitivity changes more moderately (15-30%) for medium and high 
WIMP masses due to the flatter expected recoil spectra. 

Discrimination level The presence of ER background requires a sizeable ER rejection 
factor, see also the discussion in Section 6. In the third row of Figure 6, we show how the sen¬ 
sitivity (exposure 2001 X y) depends on the discrimination level, which is varied from 99.75% 
(2.5 x 10~ 3 ER acceptance) to 99.995% (5 X 10 -5 ). For all WIMP masses, it increases lin¬ 
early with the increasing rejection level. At 99.98% rejection, it is only 20% higher compared 
to the CNNS-only case (30% NR acceptance, open symbols) for medium and high WIMP 
masses. At low masses, however, due to the dominating CNNS background, the slope of the 
improvement is weaker and the sensitivity flattens already at rejection levels >99.9%. The 
fluctuations in the curves are statistical in nature. 

Energy scale As already discussed and illustrated in Figure 3 (left), the choice of energy 
scale is important in order to achieve the best possible sensitivity. This is quantified in the last 
row of Figure 6, comparing four different energy scales: two employ a combined energy scale 
(CES), using the light (SI) and the charge information (S2) simultaneously, and two rely on 
the SI signal only, see also Section 3. Two different light yields at 122keVee, L y = 8PE/keVee 
and L y = 12PE/keVee, are used. The better resolution of the combined energy scale leads 
to significantly increased sensitivities at low WIMP masses, but also improves the sensitivity 
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for medium and - to a much lesser extent - high masses. It is interesting to note that the 
50% increase in L y has an almost negligible impact at fixed ER rejection/NR acceptance. 
Only at 10GeV/c 2 when using the Sl-only scale, the improved resolution due to the higher 
L y helps to deal with the CNNS background. In general, the increased light yield is expected 
to improve the ER rejection and therefore the sensitivity, as discussed in Section 6 and above. 

At a given Ei the Sl-only scale is somewhat more sensitive to very low WIMP masses 
(rn x < 7GeV/c 2 ) than the CES scale. This is because of its lower energy resolution, which 
allows probing the tail of the WIMP recoil spectrum by recording upward fluctuations, at 
cross sections where the CNNS background is not yet dominating. We note, however, that 
the quantification of this effect depends crucially on the relative scintillation efficiency C e ff 
employed for the study. 

A likelihood analysis was employed in order evaluate a realistic WIMP sensitivity for an 
exposure of 2001 x y and a 5-35keVnr energy interval. In energy space, ER and NR signals 
were distributed using a combined energy scale with L y = 8.0 PE/keVee, and in discrimination 
space realizing a 99.98% ER rejection at 30% NR acceptance. However, the rejection level was 
not fixed in the analysis. If no signal is observed, such an experiment will have the sensitivity 
to exclude spin-independent WIMP-nucleon scattering cross sections of 2.5 x 10 _49 cm 2 for 
40GeV/c 2 WIMPs (3.1 x 10~ 48 cm 2 at 1000GeV/c 2 ), as shown in Figure 7 (left). It improves 
upon the expected sensitivities of the XENONnT and 1.7 projects (exposure ^—20 t x y) by 
an order of magnitude, covering most of the experimentally accessible parameter space. It 
approaches the “discovery limit” suggested in Ref. [55], which is reached with an increased 
exposure of 5001 X y. 

8 Summary and Discussion 

We have studied the sensitivity of a DARWIN-like multi-ton scale LXe detector to spin- 
independent WIMP-nucleon interactions by means of toy Monte Carlo experiments, which 
take into account all known backgrounds and a realistic detector response. The simulated out¬ 
comes have been analyzed using a WIMP search box and Yellin’s maximum gap method [34], 
and by means of an unconstrained likelihood analysis for selected cases. We have examined 
the sensitivity dependence on several experimental parameters and can draw the following 
conclusions in order to optimize the dark matter reach: 

• A minimal exposure of 1001 x y is required to probe cross sections of a few xl0~ 49 cm 2 
at ~40GeV/c 2 , assuming realistic detector parameters. Sensitivities as low as 2.5 x 
10 _49 cm 2 can be reached with 2001 X y (see Figure 7, left). 

• An energy threshold of 5keVnr or below must be achieved in order keep sensitivity to 
WIMP masses below ~10GeV/c 2 . However, the CNNS-induced background increases 
significantly at lower energies. 

• The achievable sensitivity crucially depends on the rejection efficiency for ER back¬ 
grounds, which mainly stem from low energy solar neutrinos. The goal is 99.98% re¬ 
jection or better. Such a discrimination level has already been demonstrated by LXe 
experiments [8, 62]. The mean expected sensitivity at this level, taking into account all 
backgrounds, is only ~20% worse than the (unrealistic) case in which all sources besides 
the one from CNNS can be rejected completely. This factor is only mildly dependent 
on the assumed exposure and increases from 10% to 30% from 1001 X y to 12001 X y. 
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Figure 7. (left) Estimated sensitivity reach of a multi-ton LXe detector to spin-independent WIMP- 
nucleon scattering cross sections for exposures of 2001 x y (black, with lcr/2cr intervals) and 5001 x y 
(blue). It was derived using a likelihood analysis assuming a 5-35keVnr energy interval, a separa¬ 
tion of the ER and NR distributions corresponding to 99.98% ER rejection at 30% NR acceptance, 
and a combined energy scale with L y = 8.0PE/keVee. The results are compared to published lim¬ 
its from XENON100 [9] and LUX [10], and to the reach of the upcoming projects XENON1T [13], 
XENONnT [14] and LZ [15]. The “WIMP discovery limit” of Ref. [55] is also shown (red dashed), 
(right) Estimate of the typical linear dimension (diameter and height) of a cylindrical LXe time pro¬ 
jection chamber vs. the target mass used for the WIMP search. The reduction of external backgrounds 
requires a larger detector scale (we assume extra 30cm) and more LXe (right scale). 


• The concentration of nat Kr has to be reduced to a level around 0.1 ppt to achieve a 
background level which is a factor ~2.5 below the one from solar neutrinos [19]. The 
concentration of 222 Rn has to be <0.3/rBq per kg of LXe target in order to contribute 
to the background at the same level as 85 Kr. With these assumptions, the background 
is dominated by NR events from CNNS interactions at rejection levels around 99.98%. 

• The NR energy scale should be reconstructed using the light and charge signals. The 
improved energy resolution of such a scale is crucial in dealing with the CNNS events 
induced by 8 B neutrinos. 

• A high light yield is necessary to establish the required ER rejection levels, and we use 
L y = 8.0PE/keVee for most of the study. However, a higher L y does not significantly 
improve the sensitivity reach at a given rejection/acceptance level. Light yields of this 
magnitude have been achieved already, also in large detectors, e.g., LUX [61]. 

The experiment must be operated over a time span of a few years, therefore defining 
the mass-scale of the future instrument. A fiducial target of 301 of LXe would require a total 
of 6.7 y of science data to collect 2001 x y, plus additional time for detector calibration. As 
illustrated in Figure 7 (right), such a LXe target has a typical dimension of 235 cm, assuming 
a cylindrical geometry with equal diameter and height. In order to realize the required self¬ 
shielding to achieve the external background level discussed in Section 5, an additional LXe 
layer of ~15cm has to be added all around the target. This requires a total LXe mass of 431 
and an inner cryostat vessel with a typical dimension of ~265 cm. LXe TPCs of this size have 
not yet been realized, and further R&D will be necessary in order to overcome experimental 
challenges such as ultra-high voltage, target purity, high light and charge yields, charge drift, 
calibration and low backgrounds, especially in terms of 85 Kr and 222 Rn. 
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Figure 8. Sensitivity of a DARWIN-type LXe 
detector to spin-dependent WIMP-nucleon cross 
sections, assuming neutron-only couplings. Re¬ 
sults from a likelihood analysis for 2001 x y 
and 5001 x y exposures are shown (assumptions 
as in Figure 7, left), together with the limit 
by XENONIOO [67], an interpretation of the 
LUX data [ 68 ] and estimated sensitivities for 
XENON1T, XENONnT, and LZ. DARWIN and 
the high-luminosity LHC will cover common pa¬ 
rameter space. The 14TeV LHC limits for the 
coupling constants g x =g q = 0.25, 0.5, 1.0, 1.45 
(bottom to top) are taken from [ 66 ]. 


With its excellent sensitivity, a DARWIN-type multi-ton LXe detector will be a valuable 
tool to further study the WIMP particle, should it be detected in the upcoming generation 
of experiments [63]. A m x = 500GeV/c 2 (1000GeV/c 2 ) WIMP with a cross section around 
2 x 10 -4 ' cm 2 , as predicted by Ref. [64] for wino-type dark matter and at the sensitivity limit 
of XENON1T [13], would produce 46 (24) signal events in a 5-35keVnr, 30% NR acceptance 
WIMP search window of a DARWIN detector with the sensitivity shown in Figure 7 (left, 
200 1 x y). Below m x < 200 GeV/c 2 , LXe detectors are particularly well-suited to reconstruct 
the main WIMP parameters, mass and cross section, from the observed recoil spectrum [65]. 
For m x = 30GeV/c 2 and m x = 100GeV/c 2 , with cross sections of 2 x 10~ 48 cm 2 , DARWIN 
would detect 22 and 19 events from spin-independent interactions in the WIMP search region, 
respectively. We note that for WIMP masses above rs_/ 200-250 GeV/c 2 , the WIMP mass cannot 
be fully reconstructed by any target material, due to the degeneracy in po a m^ 1 for massive 
WIMPs with m x 2> mx e [63]. 

Probing identical or largely overlapping regions of parameter space is mandatory in order 
to exploit the complementarity of the different approaches to detect dark matter, with the 
goal of identifying the nature of the WIMP and studying its quantum numbers. Especially 
for spin-dependent interactions, a DARWIN-type LXe detector will largely cover parameter 
space which is also probed by the high-luminosity LHC, operating at 14 TeV center-of-mass 
energy. This is shown in Figure 8 , where we compare the LHC mono-jet sensitivity to spin- 
dependent interactions (neutron-only couplings) [ 66 ] with two multi-ton LXe exposures of 
200txy and 500txy (5-35keVnr, CES with L y = 8.0PE/keVee, likelihood sensitivity for 
99.98% ER rejection at 30% NR acceptance). The LHC 90% CL exclusion limits, reaching up 
to m x ~ 1 TeV/c 2 , were calculated for four cases, each time assuming the same values for the 
coupling constants to quarks g q and the Dirac fermion WIMP g x . The expected sensitivities 
for a DARWIN detector as well as for XENON1T, XENONnT and LZ were directly derived 
from the spin-independent sensitivities (Tsi( m \) shown in Figure 7 (left). A scaling factor 
f(m x ) was obtained by comparing the spin-independent [9] and the spin-dependent results [67] 
from XENONIOO, which are based on exactly the same dataset. The spin-dependent limits 
are then given by csd^x) = /(m x ) x Sl(, m x)- Lor LUX, we obtain the identical exclusion 
curve as published by [ 68 ] for rn x > 10 GeV/c 2 (and somewhat weaker limits at lower masses), 
confirming the validity of our approach. 

In general, the LHC will only be able to probe WIMPs up to ~1 TeV even at 14 TeV. 
For spin-independent couplings, shown in Figure 7 (left), the science reach of direct detection 
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detectors is many orders of magnitude superior to future LHC searches, which will not reach 
below ICC 44 cm 2 [ 66 ]. Therefore, direct searches - such as multi-ton LXe detectors covering 
the entire mass-space above ~6 GeV/c 2 - are necessary even at low WIMP masses to detect 
dark matter if it couples to matter via spin-independent interactions. 
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